Mon. Not. R. Astron. Soc. 000, ??-?? (20??) 



Printed 12 January 2013 



(MN IATbX style file v2.2) 



The Halo Occupation Distribution of Black Holes: 
Dependence on Mass 



Colin Degraf 1 , Matthew Oborski 1 , Tiziana Di Matteo 1 , Suchetana Chatterjee 2 , 
Daisuke Nagai 2 ' 3 , Zheng Zheng 3 , Jonathan Richardson 2 

1 McWilliams Center for Cosmology, Carnegie Mellon University, 5000 Forbes Avenue, Pittsburgh, PA 15213, USA 

2 Department of Astronomy, Yale University, New Haven, CT 06520 USA 

3 Department of Physics, Yale University, New Haven, CT 06520 USA 



o 



o 

u 

6 



> 

O 



X 



Accepted 20?? ???? ??. Received 20?? 



?; in original form 20?? xx 



ABSTRACT 

We investigate the halo occupation distribution (HOD) of black holes within a hy- 
drodynamic cosmological simulation that directly follows black hole growth. Similar 
to the HOD of galaxies/subhalos, we find that the black hole occupation number can 
be described by the form TVbh oc 1 + (AfHost)" where a evolves mildly with redshift 
indicating that a given mass halo (Mnost) at low redshift tends to host fewer BHs 
than at high redshift (as expected as a result of galaxy and BH mergers). We further 
divide the occupation number into contributions from black holes residing in central 
and satellite galaxies within a halo. The distribution of Mbh within halos tends to 
consist of a single massive BH (distributed about a peak mass strongly correlated with 
Mnost)) an d a collection of relatively low-mass secondary BHs, with weaker correlation 
with Afjjost- We also examine the spatial distribution of BHs within their host halos, 
and find they typically follow a power-law radial distribution (i.e. much more centrally 
concentrated than the subhalo distribution). Finally, we characterize the host mass 
for which BH growth is feedback dominated (e.g. star formation quenched). We show 
that halos with Mnost > 3 x 10 12 M Q have primary BHs that are feedback dominated 
by z ~ 3 with lower mass halos becoming increasingly more affected at lower redshift. 



1 INTRODUCTION 



Supermassive black holes have been found to be at th e 
center of most galaxies (|Kormendv fc Richstonej 1 1995T ) . 
and the correlation between these central black holes and 



(see, e.g. Magorrian et al.l 


19981; Ferrarese & Merrittl 


200d; Gebhardt et all I200C 


; iTremaine et al.l I2002L 


Graham & Driver 2007). Gas accretion onto supermassive 



black holes leads to AGN and quasar activity, producing 
bright objects idea for observations. One of the primary 
means of observationally studying quasars is by looking at 
their clustering proper ti es, including redshift e volution (e.g 



La Franca et al. 



r roper ti es, inc luding redshift cvofutio n (e.g 
19981; IPorciani et all 12004; ICroom et al 



20051; IShen et al.ll2007l; IMvers et al.ll2007al; [da Angel a et al 



20081 ; IShen et al.ll2009l;lRoss et al .120091 ) and luminosity de- 
pendence (see, e.g., ICroom et al.l 120051; IMvers et al.l l2007al : 
Ida Angela et al.ll2008l ; fshen et alj|2009h . By looking at the 
clustering strength of different quasar popul ations, one 
can estimate the mas s of the typical host halo (|Lidz et al.l 
l200fj ; iRoss et aUboOSl : iBonoli et ai]l2009l ; IShen et al.ll2009l ). 
thereby getting a sense of how black holes populate halos. 
This information can then be used to infe r details about 
several properties, such as quas ar lifetimes ( Haiman fc Hul 
l200ll ; iMartini fc Weinbergll200ll ), which can be constrained 



by comparing the observed quasar number density to the 
predicted density of the halo mass-function prediction using 
the typical host halos. 

Recently, quasar clustering studies have found ev- 
idenc e for a bias in th e small-scale corre l ation func- 
tion (IHennawi et al.l l200fj ; IMvers et all l2007bl . 120081 ). and 
IDeGraf et al. (20101 ) used cosmological hydrodynamic sim- 
ulations to investigate BH clustering, finding that the exis- 
tence of multiple BHs within individual galaxies has a sig- 
nificant effect on the small-scale correlation function which 
could explain this observed small-scale excess. However, this 
was an indirect means of exploring the relation between BHs 
and their typical host halos, and a more direct investigation 
is necessary to fully understand how BHs populate halos. 

With the aid of simulations, the opposite approach can 
be taken for such a direct investigation: instead of using 
clustering to predict the halos occupied by BHs, the DM ha- 
los can be directly probed to get BH occupation properties, 
which can then be extended to explore clustering properties. 
This technique has been used extensively for galaxies in the 
form of the Halo Occupation Distribution (HOD) (see, e.g. 



Berlind fc Weinbergfeool iKravtsov et~ai]|2004l ; IZheng et all 
20051) . At the most basic level, an HOD model character- 
izes the number of objects within halos as a function of the 
halo mass, and how they are spatially distributed within the 
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halo. Given these simple statistical distributions, the HOD 
mod el can be used to po p ulate halos in N-body simulations 
(e.g. iBenson et all l2000l : iBerlind et all 120031 : iBrown et all 
2008) and, assuming the occupation distribution is indepen- 
dent of the large scale en vironment (|Lemson fc Kauffmannl 
Il999l : IBerlind et al"1l2003l ) , can be used to analytically calcu- 
late the clustering statistics for a given cosmol ogical model 
(see, e.g. lSeliakll2000l ; lBerlind fc Weinberg|2002] ). The galaxy 
HOD model can then be extended further by looking at 
how the occupation properties depend upon various galaxy 
parameters, such as galaxy luminosity, color, or morphol- 



galaxy formation and evolution (e 


g. Yoshikawa et al.l 


2001 


Berlind et al.l 2003; 


Zchavi et al. 


20051; Zheng et al. 


2007 


Reid & SDereelll200£ 


: IZehavi et al. 


|2010|). 



Despite its overall success for galaxies, the HOD tech- 
nique has not been appli ed to black holes. In t his paper we 
extend the work done in IPeGraf et ail (|2010bl ) by directly 
investigating how BHs populate dark matter halos using the 
HOD formalism. In addition to characterizing the occupa- 
tion number of BHs in DM halos, we investigate the distri- 
bution of BH masses within the halo, as well as their spatial 
distribution among the component subhalos. Additionally, 
in an upcoming paper (Chatterjee et al., in prep) we will 
further extend this model to incorporate the luminosity de- 
pendencies of the black holes. By providing these details of a 
BH HOD model from a hydrodynamic simulation, we hope 
to improve the techniques available for both semi-analytic 
BH models, and theoretical studies of BH clustering. 

In Section [2] we describe the simulation used, with par- 
ticular emphasis on how the black holes are modeled. In 
Section [3] we investigate the BH occupation number both 
for halos and for central and satellite galaxies (|3.1|l , the dis- 
tribution of BH masses as a function of host halo mass (|3.2I) . 
and the spatial distribution of the BHs within their parent 
halos (|3.3[) . Finally, in section (|3.4p we look at when the feed- 
back from the BH begins to suppress further BH growth, and 
we summarize our results in Section 4. 



2 METHOD 

2.1 Numerical simulation 

In this study, we analys e the set of simulations published in 
iDi Matteo et all lliool . Here we present a brief summary 
of the si mulation code and the method used. We refer the 
reader to IDi Matteo et al.1 (|2008l ) for all details. 

The code we use is the massively para llel cosmological 
TreePM-SPH code Gadget2 (jSpringejl2o"o5r ). with the addi- 
tion of a multi-phase modeli ng of the ISM, which allow s 
treatment of star formation (|Springel fc Hernauistl [2003^1. 
and black hole accretion and associated fee dback processes 
jSpringel et al.ll2005l ; Ibi Matteo et al.ll2005l V Detailed stud- 
ies of the prescription for accretion and associated feed- 
back from massive bl ack holes and asso c iated predictions 
have b e en presented in Siiacki et al.l (120071); Di Matteo et al.l 
ll2008h: ICoiberg fc di Matted l|2008h : ICroft et alJ l|2009h : 
IPeGraf et al l l|2010ah . Important for our discussion is that 
the model has been shown to reproduce remarkably well 
both the obs erved Mbh-t relation and total black hole mass 
density pbh l|Di Matteo et al . 2008), as well as the quasar lu- 
minosity functions and its evolution in optical, soft and hard 



X-ray band l|DeGraf et al.ll2010aT l. Thus the model, within 
its intrinsic limitations, appears to serve as a fair standard 
for representing growth, activity, and evolution of massive 
black holes in numerical simulations, at least within the con- 
text of cosmological growth of black holes and not the de- 
tailed accretion physics (the detailed treatment of which is 
completely infeasible in cosmological simul ations). We note 
also that at least two indep endent teams (|Booth fc Schavel 
l2009l : ljohansson et al.ll2008l ) now have also adopted the same 
modeling for black hole accretion, feedback and BH mergers 
in the context of hydrodynamic simulations. These indepen- 
dent works, in particu lar the cosmological simulations by 
iBooth fc Schavel l|2009l ') (part of the OWLS program) have 
fully and independ ently explored th e para meter space of the 
reference model of IDi Matteo et aD l|2008h . as well as varia- 
tions to some prescriptions. This large body of already ex- 
isting work and investigations make this particular model 
somewhat of a good choice for further study. 

Within the simulation, black holes are simulated with 
collisionless particles that are created in newly emerging and 
resolved groups/galaxies. To find these groups, a friends- 
of-friends group finder is called at regular intervals on 
the fly (in time intervals equally spaced in log 10 (a), with 
A log a = log 1.25), finding groups based on particle separa- 
tions below a specified cutoff. Each group with a mass above 
5 x 10 10 /i _1 Mq that does not already contain a black hole 
is provided with one by converting its densest particle into a 
sink particle with a seed mass of Mbh, seed = 5 x 10 h~ M©. 
This seeding prescription was selected to reasonably match 
the expected formation of supermassive black holes either by 
colla pse of a supermassive star to a BH with Mb h ~ M see d 
(e.g. lBromm fc Loebll2003l ; iBegelman et al.ll2006h or by Pop 
III stars collapsing i nto ~ 10 2 Mp) BHs at z ~ 30 followed by 
expo nential growth l|Bromm fc La rson 2004; Yoshid a et al.l 
2006), reaching ~ M see( j by the time the group reaches 
~ 10 10 AfQ. After insertion, the black hole particle grows 
in mass via both accretion of surrounding gas and by merg- 
ing with other black holes. The gas accretion is modeled 

according to Mbh = ^2^,^3/2 (|Hovle fc Lvttletonl Il939l ; 
iBondi fc Ho^ll944lBondiill952t ). where p is the local gas 
density, c s is the local sound speed, and v is the velocity of 
the BH relative to the surrounding gas. Note that to allow 
for the initial rapid BH growth necessary to produce super- 
massive BHs of ~ 10 9 Mq at early time (z ~ 6), we allow for 
mildly super-Eddington accret i on (consistent with m odels 
of, e.g.. IVolonteri fc Rees 2006; B egelman et al.ll2006h . but 
limit super-Eddington accreion to a maximum of 3 x MEdd 
to prevent artificially high values. 

The accretion rate of each black hole is used 
to compute the bolome tric luminosity, L = 7?A/bhc 2 
i|Shakura fc Sunvaevll973h . Here r\ is the radiative efficiency, 
and it is fixed at 0.1 throughout the simulation and this anal- 
ysis. Some coupling between the liberated luminosity and 
the surrounding gas is expected, modeled in the simulation 
by isotropically depositing the 5 per cent of the luminosity 
as thermal energy to the local black hole kernal. This pa- 
rameter is fixed at 5 per cent based on earlier galaxy merger 
simulations such that the normalization of t he Mbh — a re- 
lation is reproduced l|Di Matteo et al1l2005l ). 

The other means of black hole growth is via merg- 
ers. When dark matter halos merge into a single halo, the 
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Table 1. Numerical Parameters 



Boxsize N p 


m DM 




c 


h~ 1 Mpc 




h- x M @ 


/i — *kpc 


33.75 2 x 486 3 


2.75 x 10 7 


4.24 x 10 6 


2.73 



iV p : Total number of particles 

mDM- Mass of dark matter particles 

m gas : Initial mass of gas particles 

e: Comoving gravitational softening length 



black holes typically fall toward the center of the new halo, 
eventually merging with one another. For these cosmologi- 
cal volumes, it is not possible to directly calculate the de- 
tails of the infalling BHs at the smallest scales, so a sub- 
resolution merger prescription is used. Since the merging 
BHs are typically found in a gaseous environment at the 
center of a ga laxy, we assume that the final coalescence 
will be rapid jMakino fc Funatol |2004| ; lEscala et all |2004 
iMaver et al1l2007T l. Thus our BHs merge when they come 
within the spatial resolution of the simulation. However, to 
prevent merging of BHs which are rapidly passing one an- 
other, mergers are permitted only if the velocity of the BHs 
relative to one another is small (comparable to the local 
sound speed). 

In addition to having been used in galaxy merger 
simulations to investigate the regulatio n of BH growth 
and correlation with host ga laxies (|Pi Matteo et ahl 



20051; 



Robertson eTall l200d : 



iHopkins et all 120073). 



Di Matteo et al. (2008) previously investigated the va- 
lidity of this method of modeling black holes in these 
cosmological simulations, finding that in addition to pro- 
ducing an Mbh — o relation that matches observations, 
the black hole mass density matches values inferred from 
the i ntegrated x-ray lum inosity function (jShankar et al.l 
2004; M arconi et al.ll20o3 ) and th e accretion rat e densit y is 
consistent with the constraints of IHopkins et all (|2007bl '). 



2.2 Simulation parameters 

The simulation analysed in this paper populates 2x486 3 par- 
ticles in a moderate volume of side length 33.75fo _1 Mpc (See 
Table [T] for additional simulation parameters). This moder- 
ate boxsize prevents the simulation from being run below 
z ~ 1 to keep the fundamental mode linear, but provides a 
large enough scale to produce statistically significant quasar 
populations. The limitation on the boxsize is necessary to 
allow for appropriate resolution to carry out the subgrid 
physics in a converged regime (for further details on the 
simulation methods, pa rameters and convergence studies see 
|Pi Matteo et al l < |2008t )1. 



2.3 Subgroup finder algorithm 

In addition to the on-the-fly friends-of-friends algorithm 
used to identify groups, a modified version of the SUB- 
FIND algorithm jSpringel et al.ll200 ll ) was run on the FoF- 
identified groups to determine the component subgroups 
{i.e. galaxies) within each group. These subgroups are de- 
fined as locally overdense, self-bound particle groups. To 



identify these regions, the algorithm analyzes each parti- 
cle within the parent group in order of decreasing density. 
For each particle i, the density of the 32 nearest neighbors 
are checked. If none are denser than particle i, it forms the 
basis for a new subgroup. If a single particle denser than 
i is found, or if the closest two denser particles belong to 
the same subgroup, particle i is assumed to be a member 
of that subgroup. If the two nearest particles denser than 
i are members of different subgroups, these two subgroups 
are stored as subgroup candidates, and are then joined into 
a new subgroup also containing i After checking each par- 
ticle in this manner, particles are checked for gravitational 
binding within their parent subgroup based on their posi- 
tion relative to the position of the most bound particle and 
the velocity relative to the mean velocity of particles in the 
group. Any particle with positive total energy is considered 
unbound, and is removed from the subgroup, leaving the 
group divided up into its component subgroups (galaxies). 



3 RESULTS 

3.1 Black Hole Occupation Number 

The most basic component of the HOD model is the oc- 
cupation number. In Figure [1] we show the mean occupa- 
tion number for both BHs (solid black line) and subgroups 
(dashed black line) as a function of host halo mass, as well 
as the the exact number of BHs found in each individual 
group (green dots). Note that these numbers are based on 
the full BH population with no mass cut; see Section 13.21 
for the BH mass distribution. We also show the contribu- 
tions to (7Vbh) arising from BHs found in the central (i.e. 
most massive) galaxy (red line) and those found in satel- 
lite galaxies (blue line). We note that this is fundamentally 
different from the traditional galaxy HOD model, in which 
'central' galaxy is of course just one. Multiple black holes 
can be found in a central galaxy (at least at low redshift 
as remnants of previous mergers) and therefore we do not 
have this restriction. For clarity, a schematic representation 
of these components (and some further subdivision we will 
discuss below) is shown in Figure [21 

In analogy with standard HOD models for the galaxy 
population we model the total, central and satellite BH oc- 
cupation number (where (iVBH,tot) = (iVBH.cen) + (-ZVBH.sat)) 
as: 



(i) 

(2) 
(3) 



respectively, where MHost is the halo mass of the host, 
Mo, Mi and M2 are normalization constants which represent 
the host masses for which we have a total of typically two 
black holes per host, two black holes in the central galaxy 
and one in a satellite galaxy, respectively. Finally a to t , ct ccn 
and a sa t are the exponents of the power law functions above. 
Note that Equations [TJ3] are not self consistent, but rather 
Equations 12131 provide an alternative parameterization from 
Equation [T] 
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Figure 1. Mean occupation number of subhalos (dashed black line), total black holes (solid black line), black holes in the central galaxy 
(red), and black holes in the satellite galaxies (blue). We also show the number of BHs in each individual group (green dots), and the 
fits to the BH occupation number using Eqns, [Tl3l and Table [2] (dotted lines). 



Standard galaxy HOD usually have a form (JVgai.tot) = 
(JVgai.cen) + (A^gai.sat } , here for BHs (JVgai.con) has been re- 
placed with a constant (equal to one) to represent our seed- 
ing condition (which artificially imposes this condition) and 
the power law forms (JVgai.sat) is similar to the ones above. 
Note again, here we have introduced an additional power law 
for modelling the A^BH.cen which allow us to characterize the 
BH numbers in central galaxies. 

Throughout this paper we refer to the most massive BH 
in a given group as the 'primary BH' while the remaining 
BHs are referred to as 'secondary BHs' (gained by merging 
with other BH- hosting halos; see Fig[3]). 

Note this makes (Af B H, SC conda ry ) = ( ) by defi- 
nition (given Equation [T]). This formulation is advantageous 
as it can be used in clustering calculati ons in the same man- 
ner as the general galaxy HOD (see, e.g. Berlind fc Weinberg! 



|2002| ; iKravtsov et aHl2004 IZhendl2004i r 

The function for (TVbh.ccii) also includes a constant 
(equal to one) since BHs are seeded in the central subgroup. 
The form of (JVBH.sat) lacks this constant since we do not 
seed subgroups with BHs, and thus there need not be any 
BHs found in the satellite galaxies. We examine each (N) to 
its appropriate form (Eqn. [Tl3)) based on halos at least twice 
the threshold mass for seeding BHs (to avoid considering 
just-seeded halos). The results of these models are plotted 
on Figure [T] as dotted lines. We emphasize that these simple 
fits are intended to provide a framework within which to get 
the typical number of BHs (total, satellite, and central) for 
the mass ranges probed in our simulation, but care should 




# Primary BH 
• +• Secondary BH 
•+«Central BH 
l •Satellite BH 

: 3 Subgroup 



Figure 2. Cartoon representation of terminology used. The grey 
circle represents a single large halo, with subhalos shown as black 
ellipses. Primary BH: the most massive BH found in the halo 
(black). Central BH: All BHs found in the central (most massive) 
subhalo (red+black). Satellite BH: All BHs found in satellite sub- 
halos (blue). Secondary BH: Any non-primary BH (blue+red). 



be taken when extrapolating to higher masses, particularly 
at high redshift where we have few data points. 

We see that (A^bh) exhibits a general trend of mildly 
decreasing slope (a) with decreasing redshift, and that ha- 
los at low redshift tend to have fewer BHs than halos with 
comparable mass at high redshift. To further understand 
this effect, in Figure [3] we show the ratio of (ATBH.secondary) 
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Table 2. Best fitting HOD parameters for Equations IU3I 



Redshift 





1 


2 


3 


5 




82 


0.90 


0.98 


1 3 


M 


1.7 x 10 12 


8.1 x 10 11 


7.2 x 10 11 


7.3 x 10 11 




1.02 


1.1 


0.94 


2.0 


Mi 


3.7 x 10 12 


4.4 x 10 12 


6.2 x 10 12 


2.9 x 10 12 


«sat 


0.49 


0.74 


0.85 


1.1 


M 2 


6.4 x 10 12 


9.5 x 10 11 


7.6 x 10 11 


7.9 x 10 11 



A 




14 



M Host ( M ©) 

Figure 3. Ratio of (NBH.sccondary) to (N snbha x ) at z=l (black), 
2 (red), 3 (green), and 5 (blue). 

to (iVsubhaios)- We note that {N suhixalos ) is sensitive to the 
mass threshold used to define a subhalo. In this paper, we 
consider any subgroup found by the SUBFIND algorithm 
(see section 2.3) to be a subhalo. However, we note that al- 
though the normalization of (iV su bhaios} is sensitive to the 
mass threshold, the slope is not. To avoid delving into the 
issue of subgroup definitions and the model-dependencies 
therein, we limit ourselves to investigating how the ratio 
evolves, which does not exhibit significant dependence on 
the mass threshold. In low mass halos, the ratio does not 
evolve with redshift, and at each redshift there are fewer 
BHs per subhalo for the low-mass halos. This is expected 
since these are halos close to the threshold mass for seeding 
a BH, and few will have undergone a merger that can lead 
to a secondary BH. However, we note that more massive ha- 
los at low redshift tend to have fewer BHs per subhalo (and 
fewer BHs in general) than the corresponding halos at high 
redshift. This decrease is likely a result of changes in the halo 
and BH merger rates. Because the number of BHs in a halo 
depends both on the rate at which new BHs enter the halo 
(via halo mergers) and the rate at which BHs within the halo 
merge with each other, a decrease in the halo merger rate 
relative to the BH merger rate would explain the decreased 
number of BHs. 

To characterize the scatter in Abh, in Figure [4] we plot 
the probability distribution P(Abh, secondary | M) (the proba- 
bility of a halo of mass M hosting Abh, secondary BHs in addi- 
tion to the initially-seeded, primary BH) for several ranges 
of host halo masses (filled circles with Poisson error bars). 
For comparison, we show a Poisson distribution with the 
same mean as the specified host mass range (dotted line). 
We see that the distribution of Abh, sat is extremely close to 
the Poisson distribution for the lower mass ranges (where 
our simulation has a large sample of halos). Even for high 



mass halos, where our statistics are poor, it appears largely 
consistent with a Poisson distribution. We note that we have 
chosen to plot the probability distribution for the secondary 
BHs rather than the total number of BHs since the existence 
of the primary BH is a condition enforced by our simulation 
which distorts P(N\M) away from a Poisson distribution 
(by removing the possibility of Abh = 0). However, this is 
an expected effect from the model for seeding BHs within 
our simulation, and we emphasize that the the physically- 
significant A^bh. secondary is well fit by a Poisson distribution 
about (Abh, secondary}- We also note that although we have 
only plotted the distribution of A^bh, secondary, we have also 
found that AaH.ccn — 1 (to avoid inclusion of the model- 
imposed primary BH) and A^BH.sat both follow an approxi- 
mate Poisson distribution as well. 

3.2 Black Hole Conditional Mass Function 

Another important facet of the model is the mass distribu- 
tion of BHs populating halos. To investigate this we pro- 
duce a conditional mass function (CMF), similar to the 
conditional lumi nosity function done for galaxies (see, e.g. 
I Yang et alfeuOSD . For this, we use -r-, — ™ — , the mean num- 
ber of BHs per logarithmic BH mass bin found in halos of 
mass MHost • We plot this quantity in Figure [5] for several 
host halo mass ranges at z=l,3. Here we see that in low 
mass halos, regardless of redshift, the BHs tend to be close 
to the seed mass, as expected. In higher mass halos, we find 
the distribution to be bimodal. The primary BH mass is dis- 
tributed about a MHost-dependent peak, while the remaining 
BHs follow an approximate power law, with the majority of 
secondary BHs being near the seed mass. The secondary BH 
CMF is also mass-dependent, with the CMF in lower mass 
halos dropping faster with increasing Mbh than in lower 
mass halos. The secondary mass distribution always peaks 
at the seed mass, however, regardless of redshift or host halo 
mass. 

The conditional mass function for the primary (i.e. 
most-massive) black hole can be approximated by a Gaus- 
sian distribution (in log 10 (A4BH)) with an integrated area of 
one 

rfA^BH.primary = 1 ^ QSSM^^tl 

d log Mbh % /2^ 2 

The best fitting parameters for this function are provided in 
Table [3] Note there is only a single group above 1O 13 M0 at 
z=3, so no fitting parameters are given for that mass range. 
From these fits we can see that the typical primary BH mass 
grows roughly proportionally to the host mass range, and 
is approximately independent of redshift. Furthermore, the 
CMF for secondary BHs is reasonably well-fit by a simple 
power law 

dN-BH, secondary _ / MbH \ 

d log Mbh ~ ^MchosJ ' U 

the parameters for which are given in Table [3] for both z=l 
and z=3. Overall, the clear trend is for the secondary BHs 
in small halos to be more strongly concentrated near the 
seed mass, while larger halos are more likely to have more 
massive secondary BHs. We emphasize that these trends are 
derived from the distributions for the mass ranges probed 
within the simulation, but are not reliable if used for masses 
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3 4 5 
'BH.secondory 



12 3 4 5 
^BH.secondory 



Figure 4. The probability distribution of A^bh, secondary for several host halo mass ranges. Filled circles show the results of our simulation 
(with Poisson error bars), and the dotted line shows a Poisson distribution centered about (A?bh, secondary) for the given host halos. 



Table 3. Fit for conditional mass function (Equations [4] and 



Host Mass Range I log 10 I -M— 





11.5 - 12.5 


12 - 13 


12.5 - 13.5 


13 - 14 


Primary BHs: 










Z = l fl 


6.68 


7.48 


8.10 


8.77 


17 


.48 


.50 


.42 


.31 


z=3 fi 


6.54 


7.28 


8.30 


N/A 


a 


.38 


.75 


.52 


N/A 


Secondary BHs: 










z=l a 


-1.51 


-1.09 


-1.01 


-1.0 


foglO ( M 0,Host) 


6.12 


6.42 


6.75 


7.21 


z=3 a 


-1.79 


-1.15 


-.71 


N/A 


foglO ( M 0,Host) 


6.29 


6.74 


7.37 


N/A 



at or below the seed mass. We also note that fi is roughly 
proportional to the host mass, so the use of finite binsizes 
in log 10 (Mhost) increases a above the ideal values for fixed 
halo masses. 



3.3 Black Hole Radial Distribution 

Another aspect of the HOD model is the spatial distribution 
of black holes within their parent halos. Although we have 
already analyzed how they populate subhalos within their 
hosts, it may also be useful to understand how BHs populate 
halos for which the subhalos have not been identified. For 
this, in Figure [6] we show the radial distribution of subhalos 



(dashed lines), BHs (solid lines), and secondary BHs (dot- 
ted lines) at redshift 1 for host groups separated into three 
mass bins (black: lO n -lO 12 M , blue: 10 12 -10 13 M Q , green: 
10 13 — 10 14 Mq), expressed as both a number density in units 
of -RJoo (top), and simply as a number per radial bin (mid- 
dle). In the lower plot we show the ratio of BHs to subhalos 
as a function of radial distance from the center. We do this 
using the complete BH population (left) and using only BHs 
above 10 7 Mq (right). Although only shown for z = 1, we 
find very similar results for z — 3, 5 as well. Thus, although 
the number of BHs found in halos of a given mass changes 
with redshift (as seen in Table [5}, the manner in which they 
are distributed within the halos remains approximately the 



Figure|6]shows that for any host mass range, black holes 
are substantially more centrally-concentrated than the sub- 
halos, as pred i cted by the black hole clustering properties in 
IPeGraf et al.l (|2010bl ). In fact, they follow a fundamentally 
different profile which can be modeled by a simple power 
law rather than the more typical NFW profile (again in 
keeping with the results of IPeGraf et al.l l2010bl ) . This in- 
creased concentration is a result of mergers between BH- 
hosting subgroups, typically between the central subgroup 
and a satellite. Because the central subgroup absorbs the 
satellite subgroup, the concentration of the subgroups does 
not increase, but that of the BHs will, since the BH will 
survive for a non-negligible time before a merging by the 
primary BH (and if it does, this is set by dynamical friction 
which is solved for in the simulations). Essentially, it is the 
existence of non-primary central BHs (see Figure [2} which 
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67896789 
Log(M BH ) Log(M BH ) 

Figure 5. Distribution of BH masses in host halos of three different mass ranges at z=l and 3. Note: Only a single group above 1O 13 M0 
exists at z=3, so that distribution has not been plotted. 



produces this increased cent r al conc entration. For further 
details see also iDeGraf et alj (|2010bl V 

To provide a means for modeling the spatial distribution 
of black holes within their parent halos (when populating ha- 
los directly rather than populating the central and satellite 
galaxies), we provide a simple fit to the radial profile (for 
r < 2 x r2oo) in form of a power law 

ri20o(r) oc (-?—) , (6) 
\r 2 oo/ 

which should be normalized to the occupation number found 
with Equation [T] The values for /3 are listed in Table U We 
note that for low-mass halos (below ~ 10 1 M©), the BHs 
are distributed roughly uniformly (in logarithmic bins) out 
to ~ i?200 5 while the most massive halos tend to have more 
BHs further from the center (though the density of BHs 
nonetheless decreases rapidly with radial distance). We also 
note that when approaching the group center, the ratio of 
BHs to subhalos tends to grow faster for low-mass halos 
than for more massive ones. Although this means that BHs 
in low mass halos are typically more centrally concentrated, 
we note that this is a result of the increased importance of 
the primary BH due to the smaller BH occupation num- 



Table 4. Radial distribution parameter for functional form 



Host Mass j3 


10 11 - 10 12 M Q 


-3.14 


10 12 - 10 13 M Q 


-2.89 


10 13 - 10 14 M Q 


-2.20 



ber. However, we note that the ratio of secondary BHs to 
subhalos (dotted lines) is constant, regardless of host halo 
mass. 

We also plot the radial distribution for massive BHs 
(Mbh > 10 7 Mq, right column) in Figure [6] Here we see 
that in addition to being less common, the more massive 
BHs are more centrally concentrated than the less massive 
BHs. This is expected since the massive BHs are almost 
exclusively primary BH, and only rarely are they satelite 
BHs (see Figured]), an d thus they should be more highly 
concentrated toward the center of the group. 
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0.01 0.10 1.00 10.00 0.01 0.10 1.00 10.00 

Radial Distance from Group Center(r/r 200 ) Radial Distance from Group Center(r/r 200 ) 

Figure 6. Top: Number density of subhalos (dashed lines), BHs (solid lines) and secondary BHs (dotted lines) in units of R^qq for all 
BHs (left) and for BHs above 1O 7 M0 (right). Middle: Mean number of BHs (solid lines) and subhalos (dashed lines) per radial bin at 
z=l. Bottom: Ratio of BHs to subhalos as a function of radial distance from group center. In all plots, color represents the host group 
mass range (black - 10 11 - 1O 12 M ; blue - 10 12 - 1O 13 M ; green - 10 13 - 1O 14 M ). 




10'' 

Holo Moss [M Q ] 



Figure 7. The fraction of groups which are feedback suppressed 
d efined as Mbh above th e quoted scatter of the M — a relation 
of lDi Matteo et alj fcOQgfi for z=l,3,5. 



3.4 Black Hole Feedback Suppression 

One final aspect of our analysis we would like to character- 
ize is the mass of dark matter halos for which black hole 
feedback has been significant (and, for example, has been 



responsible for shutting down star formation in its halo.) 
In our model, and in others in the literature it has been 
shown that feedback from a central BH can create outflows 
that are able to expel a substantial amount of the gas from 
the host galaxy, thereby suppressing further growth of the 
BH, reproducin g the M — a relation , and shutting down 
star formation (|Pi Matteo et al.|[2005l; ISpringel et al.ll2005l ; 



iHopkins et ai1l2006l . l2007al ; |Pi Matteo et alJl20oj ) 



Here we fit our simulations to provide a probabilty for 
a given mass halo to have been sign ificantly affected by BH 
feedback (i.e. feedback dominated) . |Pi Matteo et al.l (|2008l ) 
showed that our model reprduced the observed M — a re- 
lation and it does so as a result of BH feedback. We use 
the M — a relation from our simulation to obtain the black 
hole mass and the respective halo mass for which BHs are 
feedback dominated. 

In Figure [7] we show the fraction of groups in each 
halo mass bin whose BH is large enough to be considered 
'feedback suppressed' (where the feedback is strong enough 
to suppress further growth), using the condition that any 
Mbh above the quote d scatter of the M — a relation of 
|Pi Matteo et al.l l|2008l ) is feedback suppressed. We note that 
the exact choice of cutoff threshold has a mild effect on the 
overall amplitude (i.e. the exact fraction of feedback sup- 
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pressed halos), but the general trends are insensitive to the 
cutoff criteria. We find that at high redshift, very few ha- 
los have sufficiently large BHs to be feedback regulated, with 
larger halos being slightly more likely to have reached it than 
smaller halos. As time passes, the BHs in the high mass ha- 
los become more likely to become feedback regulated, and 
gradually the less-massive halos begin to become suppressed 
as well. 



4 CONCLUSIONS 

• The BH occupation number is well-described by the 
functional form (Nbh) = 1 + (jf-) tot f° r directly popu- 
lating dark matter halos. Alternatively, separate occupation 
numbers can be obtained for BHs in central and satellite 
galaxies (Eqns. [2][3]) to populate subgroups in an N-body 
simulation (or galaxy HOD model). 

• In general, (Nbh) typically follows (iV su bhaio) fairly con- 
sistently, suggesting BHs populate subhalos similarly regard- 
less of host halo mass. At low redshift, however, we find 
there are fewer BHs in the hosts (both total and relative to 
Asubhaio), particularly in moderate- mass halos, presumably 
as a result of the changing merger rates of both halos and 
the BHs within halos. 

• The scatter in (Nbh) is well described by a single pri- 
mary BH and a number of secondary BHs that follow a 
Poisson distribution about the mean secondary occupation 
number (A^bh, secondary) = ( M j5° Bt ) atot ■ We also find that the 
central and satellite occupation numbers follow approximate 
Poisson distributions. 

• The conditional mass function for the primary BH 
peaks around a BH mass strongly correlated with A/Host- 
The secondary BH mass distribution is peaked at the seed 
mass, and falls off as a power law in A/bh- The power law 
is steepest for smaller host halos, such that more massive 
halos have a wider spread of BH masses, as expected. 

• The spatial distribution of black holes within halos is 
fundamentally different from that of subhalos, tending to 
follow a power law rather than an NFW profile, leading to a 
significantly stronger central concentration of BHs relative 
to both subhalos and the underlying dark matter distribu- 
tion. Th is increased concentra tion supports the predictions 
made in iDeGraf et ail (|2010bh . though more direct investi- 
gation into our HOD-predicted correlation function will be 
investigated in an upcoming paper. 

• For a given host halo mass, the spatial distribution of 
black holes does not evolve with redshift. Thus although the 
number of BHs per host halo changes with z, how they are 
distributed within these halos remains generally the same. 

• At high redshift, few BHs are sufficiently massive to 
reach the observed M — a relation. When moving to lower 
redshifts, the more massive halos are generally the first to 
reach the M— a relation, with the lower-mass halos reaching 
the relation last. This suggests that the larger halos become 
suppressed by BH feedback at early time, and only at late 
times do the smaller halos begin to experience these sup- 
pressing effects of BH feedback. 

• We have provided best fit parameters for the mean oc- 
cupation number (Nbh) as a function of host group mass, 
as well as the BH-mass and spatial distribution functions 
within these halos to provide the necessary information to 



populate dark matter halos with black holes. Alternatively, 
we have provided the mean occupation number of BHs found 
in the central and satellite galaxies ((Nbh,csu) and (A^BH,sat), 
respectively) to provide the necessary information for di- 
rectly populating subgroups with BHs. 
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